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How achromatic is the stellar scintillation on large 
telescopes? 
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ABSTRACT 

The atmospheric scintillation of stars is the main limitation of the accuracy of ground- 
based photometry of astronomical objects. This becomes particularly noticeable for a 
study of a variability with amplitudes on the order to thousandths of stellar magnitude 
or less. We examined the problem of the colour scintillation — fluctuations of difference 
between light intensities measured simultaneously in two different photometric bands. 

Relations between the colour scintillation power (index) and the atmospheric tur- 
bulence, telescope diameter and the characteristics of the photometric channels are 
derived. Asymptotic dependencies for large telescopes (1-10 m) are obtained that allow 
to predict the value of the colour scintillation for a particular telescope and detector. 
It is shown that the colour scintillation index is oc D for measurements with both 
short (milliseconds) and long (seconds) exposures. The impact of the atmospheric 
dispersion which amplifies colour scintillation away from the zenith is estimated. We 
show that colour scintillation in the long-exposure regime depends strongly on the 
wind direction in the upper atmosphere. 
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1 INTRODUCTION 

Photometric errors due to stellar scintillati on have been 
investigated in a va riety of works (see e.g. lYound 1 19691 ; 
iDravins et al.|[l997al lbl). A distinctive feature of these fluc- 
tuations is a very small angular correlation distance, on 
the order of the isoplanatic angle (typically a few arcsec- 
onds). Consequently, the use of a comparison star at larger 
angular distance only increases the scintillation errors. In 
some cases, it the scintillation, rather than photon statis- 
tics, defines the photom etric accuracy (He aslev et al.|[l996l ; 
lEverett fc Howell [200ll ). The problem of the scintillation 
noise becomes very acute at integration times on the order 
of 1 s or less and on large telescopes. 

A variety of methods were proposed for reducing 
the scintillation noi s e: use of multiple reference stars 
|Dravins et al.|[l997a ; Gilliland et al.|[T9 93). choice of opti 



tensity simultaneously in several photometric bands (chan- 
nels) can significantly reduce the influenc e of the scin- 



tillat ion. Such approach was proposed in JCacciani et al 



mal spectr al bands 



dPravins et al.lll997bT ). arrays of small 
telescopes |Dravins et al.l 1 19981 ), conjugatio n of the tele- 
scope entrance pupil to high turbulent layers jOsborn et alj 
|2010|), selection of optimal site for te lescope installation 
(|Heaslev et al.lll996l ; lKenvon et al.ll2006l ). This multitude of 
approaches demonstrates the relevance of the problem. 

In case when the colour of an astronomical objects is 
variable, the use of an instrument which measures light in- 



E-mail: victor@sai.msu.ru 



| 2005f) to study intrinsic stellar oscillations (jRandall et al 
2005). Various types of variability in close binary sys- 
tems also often depend on the wave length of the radiation 
|Giovannelli fc Sabau-Graziatill200ll ). 

It is commonly accepted that the stellar scintillation is 
achromatic in large telescope, i.e. when the diameter of the 
receiving aperture D is much larger than th e Fresnel radius 
of scintillation tf ~ 10 cm. (|Roddier] [l981). However, this 
is true only asymptotically, but in actual situations there 
always remain some dependence on the wavelength A. 

The aim of this paper is to study the correlation of the 
scintillations in two spectrally different photometric bands. 
Typically, this correlation will be considered in terms of 
power of the differential signal — differential colour scin- 
tillation. A general description of the phenomenon is given 
in Sect.d 

In Sect. O an asymptotic dependence of the power of 
colour scintillation on telescope diameter and altitude of the 
turbulence is studied. In this section, the case of short ex- 
posures is discussed, when smoothing of scintillation due 
to wind-driven turbulence motion is negligible. The case of 
significant averaging, i.e. long-exposure regime, is analysed 
in the following Section. In both regimes, the atmospheric 
dispersion greatly enhances the colour scintillation. 
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In the last Section, the effectiveness of multicolour tech- 
nique (in the sense of simultaneous light intensity measure- 
ment in different photometric bands) for the scintillation 
noise suppression is evaluated. 



2 GENERAL FORMULAE 

Basic expressions to calculate the covariance of light inten- 
sity fl uctuations in tw o photometric bands are given in the 
paper (|Kornilovll2011al ). They were obtained for the approxi- 
ma tion of weak per turbations by the method previously used 
m (|Tokovininll2003l ). The starting relation for the covariance 
01,2 of normalised intensities, arising due to turbulence in 
the layer located at the distance z from the receiver, is the 
expression: 

Cl,2 = {hh)-{h){h) « 

« if <Xi(Ai)x2(A 2 )> F 1 (A 1 )F 2 (A 2 ) d\ 1 d\ 2 , (1) 

where Xi an d X.2 are the averaged over entrance aperture 
logarithms of wave amplitudes in the first and the second 
channels, and Fi(A) and -p2(A) are normalised products of 
the photon sensitivity in the photometric bands with the 
energy distribution of the source spectrum. This equation 
is approximate because the right-hand side is written for 
the log-amplitudes and the left-hand side for the intensities. 
For weak scintillation, the expression can be regarded as 

accurate. 

Using the formulae 7-11 from the paper jTokovininl 

2003), one can calculate the covariance of the log-amplitudes 
(Xi(Ai)X2(A2)) where apertures of the photometric channels 
may be not be coincident, in general. For coincident circu- 
lar apertures, the product of the Fourier transforms of the 
aperture functions W1W2 is replaced by the axi-symmetric 
aperture filter A(f) during further integration over the two- 
dimension spatial frequency. 



2.1 Correlation of the scintillation in two 
photometric bands 

The relationship between the scintillation power (scintilla- 
tion index) s 2 and the turbulence intensity C 2 (z) at the dis- 
tance z is described with help of weighting functions (WF), 
that is scintillation power produced by a turbulent layer of 
unity intensity. In the linear approximation of the weak per- 
turbations 



s 2 = / C 2 (z)Q(z)dz. 



(2) 



where Q(z) is the respective WF and the integration is 
performed over whole atm ospheric depth. In analogy with 
the expression for Q(z) in l|Tokovininll2003T ). we can obtain 
an expression for the covariance weighting function Ri^(z) 
which describes the covariance of scintillation in two differ- 
ent photometric channels: 

/>oc 

R l<2 (z) = 9.61 / r B/3 Si,2{z,f)A(f)df, (3) 
Jo 

where the integration is performed over the modulus the spa- 
tial frequency /, A(f) is the spectral aperture filter, equal to 



A(f) = (2Ji(nDf)/(irDf)) 2 for the case of coincident cir- 
cular apertures of the diameter D, and the kernel Si, 2(2, f) 
is the Fresnel spectral filter. This filter is a double integral 
which splits into a product of two integrals: 



Fi(A 



Ai 



— sin(-7rAi2:/ 2 ) dAi 



F 2 



^^-Sm(TT\2zf 2 )d\2. (4) 

A2 



Then, the observed covariance ci,2 of the scintillation 
in two channels (cross-index) depends on the vertical distri- 
bution of the structure coefficient C„(z): 



Cl,2 



Cl(z)R 12 (z)dz. 



(5) 



If D S > rjr = (A2) 1 / 2 ), we can use the method 
of i|Roddierl Il98ll ). putting sin(7rAz/ 2 ) w TvXzf 2 . The 
asymptotic dependence for the covariance WF R\,2(z) = 
17.33 D~ 7ri z 2 is obtained. This expression is identical to 
the known relation for the scintillation power in the case of 
large telescopes. 



2.2 Colour scintillation 

Correlation of the scintillation can be described in terms of 
the differential colour index — the variance of the difference 
of normalised signals in two photo metric channels, an ana- 
logue of differential spatial index l|Tokovininl |2002| ) . In the 
case of strong correlation expected for large telescopes, this 
approach is more practical for measurement and analysis. 
This colour scintillation arises from the dependence of the 
diffraction on the wavelength of light. Based on the relation 
for the variance of a difference: 



Sd = s\ + s 2 — 2C1,2, 



(6) 



the colour weighting function Qd(z) can be represented as: 

Q d (z) = Qi(z) + Q 2 (z)-2R 1)2 (z), (7) 

Returning to the expression (|3]) and summing the corre- 
sponding integrands, we obtain 



Q d (z) = 9.61 / /- 8/s Ti, a (z, f)A(f) df 
Jo 

where Ti j2 (z,/) is the differential Fresnel filter. 



(8) 



Ti, 2 (z,f) = J J dAadAi 

Fl,F 2 

/Fi(Ai) 2 \ 
x [ j- sin( 7 rA 1 z/ / ) 



^sinM 22 / 2 )) 2 . (9) 



If the usual Fresnel filter used for the Q(z) calculation has 
behaviour oc f 4 near F = 0, the differential Fresnel filter 
T(z, f) is proportional to f 12 there and suppresses the low 
frequencies passed by the aperture filter very effectively. 

For apertures of arbitrary size and wide photometric 
bands, the functions Ri :2 (z) and Qd(z) can be obtained 
by numerical integration with formulae <(3j , (0 and ([S]!, 

as done when calculating the WFs in the program at- 
mos developed for processing of meas urements with Multi- 
Aperture Scintillation Sensor (MASS) jKornilov et al.ll2003l : 
iTokovinin et al.ll2003l ). The computed functions Qd(z) for 
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Figure 1. WFs for the 0.3 m, 0.6 m, 1.0 m and 2.5 m tele- 
scopes with central obscuration 0.3 in photometric bands B (thin 
solid line) and R (thin dashed). The colour WFs Qb,r( z ) are 
shown with a thick line. The cases of the apertures without cen- 
tral obscuration are shown by dotted lines. The thick dashed line 
represents the estimate calculated with the formula i!3l 



the B and R photometric bands (effective wavelengths \b ~ 
440 nm Xr = 690 nm) and for four different telescope aper- 
tures are shown in Fig. [T] 



3 SHORT EXPOSURES 



In the case of large apertures D 3> (Az) ' , the asymptotic 
dependencies for both index and cross-index are completely 
achromatic (|2.1[) . Therefore, it is not possible to obtain the 
asymptote of the colour index s 2 , from the formula How- 
ever, the asymptote can be found from the analysis of the 
colour WF ©. 

Writing this formula for two monochromatic channels 
Ai and A2 and a circular aperture of diameter D and passing 
to the dimensionless frequency q — fD, we get 



) d (z) = 9.61 D 



i/3 



poo 


"2Ji(7n?)" 


/<) 





3/3 



/sin(7rAizg 2 /D 2 ) sm{-K\ 2 zq 2 / D 2 ) 



Ai 



A, 



(10) 



Clearly, the integral still depends on the parameters z, 
D, Ai and A2. In this case, the aperture A(q) and the dif- 
ferential T(q) filters almost do not overlap (Fig. [2]). In the 
region where T(q) is significantly different from zero, one 




Figure 2. Position and typical form of the aperture A(q) (thin 
line) and differential Frcsncl T(q) (thick line) spectral filters. The 
approximation of the aperture filter is shown by the dashed line. 



can use the asymptotic form of the aperture filter: 



2Ji(nq) 



nq 



8 2, 3 . 
— 5— a cos (irq — -IT) 
7T 4 q 3 4 



(11) 



which is valid when g > 1. This approximation is shown 
in Fig. [2] and one can see that in the area of interest, it 
coincides with the exact function A(q). Since A(q) oscillates 
within the envelope of T(q), it is possible to replace cos 2 
with its average value 1/2: 



Q d (z) = 9.6173 



5/3 



-17/3 



sm(irXizq 2 /D 2 ) sin(7rA 2 2<? 2 /-D 2 



Ai 



A, 



dq, (12) 



This operation allows us to reduce the integrand to 
a dimensionless form by replacing the variable u = 
TviX^^zD-iq 2 : 

Q d (z) = 2.85 z 7/3 D-' A (X 1 X 2 ) 1/6 



-10/3 



— sin(em) — a sin(u/a) ) du, (13) 



s a 

l ' 2 . The integral Xi in this formula can 



where a = (A1/A2 
not be represented as a sum of three integrals, because each 
separate component diverges. We could not calculate X\ an- 
alytically, therefore its behaviour as a function of the ratio 
A1/A2 = a 2 was calculated numerically and is presented 
in Fig. [3] Near the point of Ai = A2 the dependence is 
quadratic, both branches increase indefinitely at small and 
large arguments. 

3.1 Analysis of the asymptotic dependence 



The obtained asymptotic dependence (|13[) is confirmed by 
the numerical integration of the formula {SJ. The polychro- 
matic WFs for 0.3 m, 0.6 m, 1.0 m and 2.5 m apertures and 
B and R photometric bands are shown in Fig. [T] Complex 
structure of T(q) affects the form of the WFs, they have 
a noticeable waves with altitude. The figure also shows the 
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Figure 3. Integrals X\ in the expression [13] (solid line) and X 2 in 
the expression 1241 as a function of the parameter a 2 = A1/A2. 



WFs for the same aperture without a central obscuration for 
comparison with the asymptotic estimates (|13l) . In the scale 
of the figure, the both curves coincide to within 1-3 percent. 

Such a good match between the polychromatic and the 
monochromatic curves occurs because the scintillation is 
only slightly chromatic in the case of a large aperture. As a 
consequence, the functions Qd(z) are virtually independent 
of the spectral band width. This means that estimates based 
on the formula (|13|) can be used for real photometric bands. 
To account for the central obscuration 0.3 (typical for a two- 
mirror telescope), it is sufficient to multiply the asymptotic 
estimate by the factor 1.9. For the normal scintillation, this 
factor is less than ~ 1.21. 

Evident relation between the altitude h and the distance 
z — hM z predicts that the colour scintillation power is a 
function on the air mass M z in the form of s 2 d oc Mi ^ 3 due 
to simple geometric factor of light propagation. 

The dependence of colour scintillation on altitude is 
stronger than for the conventional scintillation. It means 
that the colour scintillation is generated at the higher layers 
of the turbulent atmosphere. Given that the scintillation in- 
dex for t he case of sho rt exposure is defined by the known 
relation (|Roddierlll98lh : 



s 2 = 17.33 D~ 7/3 [ Cl(z)z 2 Az, 

J A 



(14) 



and that the colour index is equal to 



4 = 2.852 L»- 3 (A 1 A 2 ) 1/6 Xi(A 1 /A 2 ) / C 2 n {z) z 7/3 dz, (15) 

J A 

it is possible to estimate the relative power of s d / s 2 — achro- 
matic factor of the scintillation. For that, we take z 1 ^ 3 out 
the integral, replacing it with average value (z 1 ^ 3 ) weighted 
with C 2 (z) z 2 . Then the ratio of the indices can be written 
as 



% = 0.165 £>- 2/3 (AiA 2 



,V6/„V3 



}Zi(Ai/A 2 



(16) 



Estimates of the altitude atmospheric moments M p = 
f C 2 (z)z p dz can be obtained from vertical distributions 
of the optical turbulence m easured with the MASS instru- 
ment (|Kornilov et al.l 120031 ). For example, measurements 



at Mount Shadzhatma z in 2007-2009 show the median of 
M 2 = 1.1 • 10" 5 m 7/3 (|Kornilov et al.l 120101 ). The quantity 
(z 1//3 ) depends only weakly on the actual turbulence profile 
and can be set equal to 24 m 1//3 . This assumption ensures 
accuracy of the estimate not worse than 10 percent. 

The expression (|16p shows that the achromatization 
comes not very quickly: increasing the telescope diameter by 
2.5 times reduces the chromaticity only by half (the scintilla- 
tion amplitude is reduced by 1.5 times). For small apertures 
(on the order to the Fresnel radius), the colour scintillation 
is ~ 0.1 of the normal one (or ~ 0.3 in amplitude). In the 
case of 2.5 m telescope, the achromatic factor is ~ 0.003 
in power or ~ 0.05 in amplitude when measurements are 
performed in the B and R bands. 



3.2 Decorrelation due to atmospheric dispersion 

When measurements are done away from the zenith, beams 
falling into the first Ai and the second A2 channels, pass 
through atmosphere by a few different ways because of the 
effect of atmospheric dispersion. That leads to a decorrela- 
tion of the scintillation, i.e. to an increase of t he colour scin- 
tillation. This phenomenon was discussed in |Dravins et al 



1997bl : lHubbard fc Reitsemalll98ll;l Caccia. Vernin fc Azouiti 



In the paper ( Kornilovl 201 lal ). the expression for the 



additional power As 2 of the colour index s d was presented. 
The WF for calculating the extra power due to atmospheric 
dispersion is described by the formula: 

00 

Qa(z) = 9.61 j /- 8/3 S!, 2 ( 2 , /) (2 - 2J (2nxf))A(f) df, 


(17) 

where A(f) is the aperture filter as usual, and Si,2(z,/) is 
the Fresnel filter of the kind of Q, and x — x(z) is the 
displacement of one beam with respect to the other at the 
distance z. 

Recall that for a flat isothermal atmosphere, the beam 
displacement x at the altitude h is given by: 



x — tan 7 sec 7 Ani, 2 e h °l H ° (Ho — (Ho + h)e 



-h/H \ 



(18) 



where 7 is the zenith distance, h is the observatory level 
above sea, h = z cos 7 is the height of turbulent layer with 
respect to the telescope level, Ho is the atmospheric height 
scale (~ 8300 m). The multiplier can be written also as 
tan 7 sec 7 = M z (M 2 — l) 1 ^ 2 if the air mass M z is used 
instead of 7. We emphasise that the x is not a horizontal 
shift, but displacement perpendicular to the light direction. 
With height increasing, x grows approximately quadratically 
and at altitude about 14 km (1.68 Ho) it reaches 1/2 of its 
maximum value, then asymptotically approaches it. 

It is convenient to estimate the effect when tan 7 sec 7 = 
1 or air mass M z — 1.27202 = \f$ (<f> is the golden ratio). 
The shift outside the atmosphere is 1™ = Ho Aru, 2 e~ h °^ H ° . 
For the B and R photometric bands and observatory at 
2500 m, the difference between the refractive indices at the 
sea level is Aub,r = —5.15 • 10 -6 and Xoo = 0.032 m. 

Let us investigate the Qa(z) behaviour in the case of 
two monochromatic channels and circular telescope aper- 
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Figure 4. Relative displacement of the beams y = x/rp versus 
distance to the turbulent layer for different air masses (labels on 
the curves) for the B and R photometric bands. 



Figure 5. Dependence of the integral X3 on the relative displace- 
ment. The initial part of the curve is shown with the solid line 
with coordinate axes on the left and below. The total dependence 
is shown with dashed line with axes at the top and right. 



ture: 



9.61 / / 

/o 



_ 8 /3 sin(7rAi2:/ 2 ) sin(7rA 2 z/ 2 ) 
Ai A2 

2 



x(2-2J (2 TO /))( g^^ j d/. (19) 

One can introduce the dimensionless frequency q — rpf and 
the relative shift y — x/tf- Calculation with (|18[) shows that 
for all reasonable air masses the condition y < 1 holds be- 
cause x increases with altitude and t\f increases with dis- 
tance to the layer. The dependence of y on the distance is 
shown in Fig. 2] for different air masses. For B and R photo- 
metric bands, the maximum of y is related to the air mass as 
0.22 Ml /2 (M 2 - 1) 1/2 . The average compound Fresnel 



y 

radius tf 



,1/2 



(AiAj 



a/4 



is used here. The main peak of 



the integrand lies in the domain q < 1. Note that the Fresnel 
cross-filter for some frequencies may be negative. 

We replace the aperture filter A(f) with its approxi- 
mation (|11[) which is applicable when q > tf/hD. Just as 
before, we can replace cos 2 with 1/2. Unlike the case consid- 
ered in Sect. [3] the cross-filter is proportional to ~ q 4 near 
zero and it does not provide effective suppression of the ap- 
proximation error. Thus, the approximate integrand has an 
infinite derivative at q = 0, while the exact expression has 
zero derivative. Such distinction leads to ~10 percent error, 
which was estimated from a comparison of the asymptote 
with numerical integration of the formula (|19|l . 

As a result of these transformations, the expression 
reduces to the dimensionless integral T^{a,y), where a = 
(A1/A2) 1 / 2 , with the dimensions factor in front of it: 



a(z) = 0.789 D"V /3 (AiA 2 ) 1/6 

q~ 17 ' 3 sin(7raq 2 ) sin(7rg 2 /a)(l — Jo(27ryq)) dg. (20) 



The I3 depends implicitly on the distance to the layer z 
through y, but it is not very important for a general esti- 
mate. 

The dependence on y can be obtained by replacing the 



filter 1 — Jo(2ivyq) with its approximate expression ir 2 y 2 q 2 
which is applicable when y < 1. In this case, the integral 



X 3 « 17.3y 2 (A 1 A 2 )- 4/b ((A 1 + A 2 ) 4A1 - (Ax - MY 1 *). (21) 

The relative displacement y depends implicitly on the wave- 
length ratio a through the refractive index difference (|18[) . 
As the difference between the wavelengths gets larger, the 
factor depending on a in (|21[l becomes smaller, but this is 
compensated by the increased y 2 . 

In the case Ai = A2, the integral I3 is calculated 
analytically after differentiation over y, but in a rather 
complicated way through hypergeometric functions. Its be- 
haviour is shown in Fig. [5] in two scales. We see that for 
y < 0.3 the dependence is quadratic, then it progressively 
becomes weaker than linear. Comparison with the formula 
(|21|l demonstrates that the latter provides estimates with an 
accuracy of about 5 percent up to y = 0.3. Consequently, 
for air masses M z < 1.5, the use of quadratic dependence of 
As 2 oc y 2 or As 2 oc M z (M 2 — 1) is sufficient for practical 
purposes. 



3.3 Relative contribution of the additional colour 
scintillation 

Comparison of the formula (|20|) with the (|13[) for the colour 
scintillation shows that their ratio depends only on the wave- 
lengths of the photometric channels and on the relative dis- 
placement y. The relative contribution of the atmospheric 
dispersion to colour scintillation will be determined as fol- 
lows: 



A^ 

o2 



= 0.277 



Vl/.Ai/Aa)) 



(22) 



Xi(Ai/A 2 ) ' 

Using the plots in Fig. [3] and Fig. [4] one can estimate that 
the ratio Asj^/s 2 , reaches 1 at y ~ 0.1. Such value corre- 
sponds to observations at air mass M z ~ 1.1. Consequently, 
colour scintillation caused by the propagation will dominate 
over colour scintillation caused by the atmospheric disper- 
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z, km 

Figure 6. Colour WFs for 1 m and 2.5 m telescopes. The func- 
tions Qd(z) + Qa(z) for short exposures at air mass 1.0 (thick 
lines), 1.1, 1.2 and 1.3 (bottom to top). 



sion only near the zenith (7 < 25°), when measured in the 
B and R bands. 

We improve the expression ()16[) by adding extra colour 
scintillation induced by the atmospheric dispersion: 

+ „ AS ^ « (0.005 + 0.021 M z (M 2 - 1)) D~ 2/3 . (23) 

Numerical integration of the initial expression taking in ac- 
count all altitudinal dependencies, leads to the WFs which 
are shown in FigdJ] where a family of functions Qd(z) + 
Qa(z) is plotted for 1 m and 2.5 m telescopes and air masses 
1.0, 1.1, 1.2 and 1.3. 



4 LONG EXPOSURES 

In the previous sections, the case of measurements with tem- 
poral scales comparable to the atmospheric coherence time 
was considered. Conventional photometry rarely has to deal 
with such short exposures and therefore, the study of the 
colour scintillation on time scales of seconds and more has 
a larger practical interest. As a rul e, on such time scales 
the long exposure regime is realised l|Kornilovll2011bl ) when 
the turbulence drifts during exposure time r at the distance 
wt 3> D (here w is the wind speed). The scintillation power 
determines the photo metric errors in this case. 

It was shown in (|Kornilovll20"Tlal ) for the case of an in- 
finitely small aperture that the wind averaging increases the 
correlation of scintillations, i.e. reduces the colour scintil- 
lation. Naturally, in that case the condition wt 3> rp was 
used as a criterion of the long-exposure regime. 



4.1 Wind smoothing 

Using the approach bas ed on the hypot hesis of frozen turbu- 
lence and described in (|Kornilo v 2011b), we modify suitably 
the expression (|10[) by adding the spectral filter T\ (wt/) of 
the wind shear, equal to 1/ttw(z)t/ in the long-exposure 
regime. Note that w(z) is the component of the wind vec- 
tor transverse to the line of sight. Omitting the intermediate 
transformations, we write the analogue of ()13|) for the colour 
WF: 

U' d (z) = 1.609 z 17/6 D- 3 (AiA 2 ) 5/12 

x / u~ 23,/6 I i sin(ow) — a sin(u/a) J du, (24) 
Jo \ a J 

The dimension of this function is m 2 / 3 because for calcu- 
lation of the scintillation index s 2 ,, it should be divided by 
the wind shear w(z)t. The integral I2 in this formula can 
be calculated numerically. Its behaviour depending on the 
ratio A1/A2 = a 2 , is shown in Fig. [3] The behaviour is not 
fundamentally different from that of X\ in (|13[1 . 

Attention is drawn to an almost cubic dependence on 
the altitude. This means that for long exposures a strato- 
spheric turbulence at altitudes of 25-30 km will contribute 
to the colour scintillation. Numerical integration of the ini- 
tial expression confirms the asymptote in the form of (I24| l. 

In the long-exposure case, the scintillation becomes 
more achromatic with increasing D faster than at short 
exposures. Indeed, using the asymptote fo r U'(z) = 
10.66 L>" 4/3 z 2 from the paper jKornilovll20lTbT ). the follow- 
ing expression can be written (the scintillation indices for 
long exposures are denoted as s 2 and s 2 ): 

4| = 0.151 D- 5/3 (AiA 2 ) 5/12 (2 5/6 )22(Ai/A 2 ). (25) 

The quantity (z 5 ^ 6 ) is the z 5 ^ 6 averaged with the weight 
C 2 (z) z 2 /w(z). The wind velocity dependence on altitude 
w(z) additionally increases the effective altitude of colour 
scintillation above the tropopause, because a decrease of tur- 
bulence intensity with height is offset by a decrease in the 
wind speed above the tropopause. Based on the actual ver- 
tical distributions of turbulence (see Sect. I3.1[l we can set 
(z 5 / 6 ) « 3300 m 5 / 6 . This estimate has uncertainty about 
20 percent. For 2.5 m telescope the s 2 d /s 2 is equal to 0.0001 
in power or 0.01 in amplitude, which is considerably less 
than in the case of short exposures. 

The dependence of the colour scintillation on the air 
mass M z is almost the fourth degree and when the wind 
blows along line of sight, becomes oc M 29 ^ 6 because the 
s 2 , is inversely proportional to the transversal wind which 
decreases as M" 1 . The achromatic factor (|25|l increases as 
Mz^ 6 at all times. 

4.2 Decorrelation due to atmospheric dispersion 

In the long exposure regime, additional colour scintillation 
due to atmospheric dispersion will depend on the angle 9 
between the direction of the wind vector and the azimuth 
of the star. Indeed, during long exposure the wavefront is 
averaged along the wind direction over an area of about 
wtD. When the displacement x is transverse to the wind, 
the non-overlapping area is ~ 2xwr, much larger than ~ 
2xD for the case when x is parallel to the wind. 
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Therefore, during the transition from two-dimensional 
integration over spatial frequency to the integration over 
its modulus, we must average over the angle <f) the product 
of the two filters, the one for aperture displacement and 
another for the wind shear. In the equation (JTTJ) the factor 
of 2 — 2 Jo (2-Kxf) is the result of the averaging over the 
angle of the asymmetric filter of the aperture displacement 



2 — 2 cos(2ivxf cos < 



4 sin (irxfcost 



We include the 

filter of the wind shear in the form of sine 2 (tor/ cos <j>) by 
rotating the coordinate x-axis by the angle in the direction 
of the wind. The product is averaged over <f>: 

r 2n 



C 



1 

2^ 



4sin (nxf cos(0 — 9)) sine (wt/cos< 



The asymptote of the integral C at wt 



(26) 

oo depends on 



the angle 9. Since xf = yq < 1, the linear approximation 
sin 2 (7ra;/cos(0 — 9)) « n 2 x 2 f 2 cos 2 (<j> — 9) may be used. For 
9 = (same wind and source azimuth), this gives: 



A = 



2 W 2 T 



Jo^nwrf)). 



(27) 



If 9 = 7r/2 (wind direction perpendicular to the atmospheric 
dispersion), the averaging results in: 



C± = 4n 2 x 2 f 2 Ti(wrf) 



2w 2 i 



r (l- Jo(47Iw/)) 



- 2 / 



4tt — (28) 



where 71 (wr/) is the wind filter l|Kornilovl 1201 lbl ) which 
equals l/irwrf in the long-exposure regime. 

Comparison of the relations for £± and C\i shows that 
both components are approximately equal at very low fre- 
quencies / < 1/ttwt. At larger frequencies, £± ^> C\\. 
Taking into account that the maximum of the scintillation 
is at / = 0.555/ri? ~ 5 — 10 (maximum of the product 
q~ 8 ^ 3 sin 2 (7rq 2 )) and that 1/ttwt <C 0.1, it can be argued 
that the investigated effect is mainly determined by the per- 
pendicular component of wind. 

Corresponding WF U'& (z) can be obtained by replacing 
the factor 2 — 2 Jo (2iryq) with the expression for C in the 
formula (|17|) or in formula (| ltjp in the case of two monochro- 
matic bands and a circular aperture. To get a formula like 
(I24[) . we have to isolate the factor 1/wr. 

Except for the factor 1 — Jq(Attwt f) which almost every- 
where is equal to unity, the component C\\ does not depend 
on the frequency and can be taken out of the integral. In 
this case 



2 2 

&a(*)|| = Ri,a(z) = 8.66— D~ 

IWT WT 



7/3 z 2 



= 8.66^-D- 7/3 z 3 (X 1 \ 2 ) 1/2 , (29) 

WT 

where Rifi(z) is the cross-index WF ([3]) having the usual 
asymptote. 

In calculating the component C±, one can neglect the 
term £« . After replacing / with the dimensionless frequency 
q, replacing the aperture filter with its approximation (|lip . 
separation of 1/wr and easy transformations we obtain: 

U' A (z)± =4.96y 2 £>- 3 (AiA 2 ) 5/1 V 7/6 

^-14/3 gijj^Qiq, 2 ) s in(7rq 2 /a) dq. (30) 




Figure 7. Colour WFs for 1 m and 2.5 m telescopes. The func- 
tions U' d (z) + U'^(z) for measurements with long exposures; the 
thick line corresponds to observations at zenith, thin solid lines 
— at the air masses 1.05, 1.1, 1.2 and 1.3 for the transversal wind 
and the thin dashed lines — at air masses 1.3 and 1.5 for the 
longitudinal wind. 



The approximation of the integrand near zero frequency is 
poor because the power of the q~ 14 ^ 3 term is close to the 
limit for integrability. Numerical integration of the original 
formula shows that the error of the asymptotic expression 
(|30[1 is on the order 20 percent for D = 2—4 m and z ~ 4 km. 
In this equation, the integral can be calculated analytically 
and it equals to: 



1' 4 = 7.175 (AiA 2 ) 



-11/12, 



((Ai+A 2 ) 11/6 -|A 1 -A 2 | 11/e ). (31) 



ll/6\ 



Relative contribution of the wind shear to 
colour scintillation 



4.3 



The formulae (|29p and (130[) allow us to calculate the relative 
increase of the colour scintillation caused by the turbulence 
wind drift. To do this, we can use the method of Sects. [3~T1 
and 14. II In the case of parallel wind we have: 



^ = 5.38 V 



(AlA; 



a/12 



wt I 2 (Ai/A 2 



■7J 2/ V /6 >, 



(32) 



where (z ' ) is the average of z 1/b with the weight 



,17/6 



/w(z). This weight has maximum at even higher 
altitudes than in the case of (|25p : for the typical conditions 
considered here we estimate (z 1//6 ) ~ 5 m 1 ^ 6 . For the B and 
R bands we can use the approximation « 30 y 2 D 2 ^ 3 /wt 
which gives the magnification factor < 2 for y = 0.3 
(M z ~ 1.5), wind shear wt = 5, and 1 m telescope. 

The proportionality to D 2 ' 3 looks strange, but one 
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should take into account that the wind shear wt must in- 
crease together with D to preserve the validity of the long- 
exposure regime. 

An entirely different situation arises in the case of per- 
pendicular wind, when the magnification of the colour scin- 
tillation due to wind depends only on the wavelengths, as in 
the case of 

s d X 2 (Ai/A 2 ) 

It is seen that the intensity of colour scintillation is tripled 
even at y = 0.05 or M z « 1.05, increasing by w 30 times at 
M z = 1.3. 

The above estimates are confirmed by the exact calcu- 
lations, the results of which are shown in Fig. The func- 
tions U'&(z) are calculated for telescopes of 1 m and 2.5 m 
diameter with central obscuration 0.3, polychromatic B and 
R bands, and the exact dependence of linear atmospheric 
dispersion x on the altitude. 

A large difference between these two situations is seen. 
The effect of the parallel wind at M z — 1.5 is much less 
than the effect of the transverse wind at M z — 1.05. Par- 
allel WFs are calculated for a small wind shear wt — 5 m; 
at larger wind shears the contribution of the atmospheric 
dispersion will only become smaller. For the perpendicu- 
lar wind, the functions do not depend on the wind shear, 
although of course, the scintillation index will decrease in 
proportion to ~ 1/wt in any case. 

In analogy with (|23p . we can write: 

As ^ ± ^ » (0.0003 + 0.012 M Z {M 2 Z - 1)) D~ 5/3 . (34) 

For this approximation, a worst-case scenario is 
adopted: the wind is perpendicular to the beam displace- 
ment due to the atmospheric dispersion. 

In the formula (|34[) as well as in the expression (|30[) . 
the wind speed is not present explicitly. It enters through 
the condition wt 3> D. This means that for any small per- 
pendicular component of the wind but rather long exposure, 
the worst case situation may be encountered. 



5 DISCUSSION 

As an example, the magnitudes of the normal and the colour 
scintillations are presented in Table. Q] for different situa- 
tions. They are calculated using the exact WFs (Fig. [6] [7J 
and a real turbulence profile from MASS m easurements at 
Mount Shatdzhatmaz (|Kornilov et al.lboioT ). One specific 
turbulence profile was chosen for which the atmospheric mo- 
ment M 2 corresponded to its median. In the long-exposure 
regime, the values s 2 wt (see Sect. 14.1] ) are given, from which 
the scintillation indices can be obtained by dividing those 
values by wt. 

These estimates confirm that colour scintillation is 
much less than normal scintillation when observing near 
the zenith. This is especially noticeable in the long-exposure 
regime typical for conventional photometry. The effect of at- 
mospheric dispersion favours observations near the zenith, 
but colour scintillation is reduced at large telescopes, allow- 
ing observations at larger air mass. 



Table 1. Estimation of the normal (N) and the colour (C) scin- 
tillation indices in the B and R spectral bands for 1 m and 2.5 m 
telescopes in the short (SE) and long (LE) exposure regimes. For 
the LE regime, two situations are considered: with the wind per- 
pendicular (LEx) and parallel (LE11) to direction of the atmo- 
spheric dispersion. For clarity, the amplitudes Vr are written in 
millimagnitudes in the SE case and in millimagnitudes per meter 
in the LE case. 



SE LEx LE|| 

M z 1 m 2.5 m 1 m 2.5 m 1 m 2.5 m 



N 


1.0 


14.48 


5.21 


10.85 


5.91 


10.85 


5.91 


C 


1.0 


1.26 


0.33 


0.23 


0.06 


0.23 




C 


1.05 






0.43 


0.12 






c 


1.1 


1.64 


0.43 


0.57 


0.16 






c 


1.2 


2.08 


0.54 


0.88 


0.25 






c 


1.3 


2.56 


0.67 


1.18 


0.33 


0.25 


0.07 


c 


1.5 










0.27 


0.08 



Obviously, the diameter of telescope aperture is also im- 
portant for another reason. The gain from potential scintilla- 
tion noise suppression by multicolour photometry is achiev- 
able only when the photon noise is low enough. Estimation 
of the limiting object is straightforward. Only two points 
will be emphasised here. 

The ratio of normal scintillation to the photon noise de- 
pend on the exposure time. In the short-exposure regime, the 
scintillation noise is oc D~ 7 ^ 6 and decreases slightly faster 
than the photon noise oc D _1 with increasing telescope di- 
ameter D. In the long exposure regime, the scintillation de- 
creases as oc D~ 2 ^ 3 , much slower than the photon noise. 
Therefore, in this situation, the scintillation noise becomes 
dominant at large telescope apertures. 

The colour scintillation, on the other hand, always falls 
faster (oc D~ 3//2 ) than the photon noise. Consequently, when 
the telescope diameter increases, the photon noise becomes, 
sooner or later, the main limitation of the measurement ac- 
curacy. So, we can take advantage of the small colour scintil- 
lation only by choosing a bright enough astronomical object. 

A similar situation exists in the c ase of photometry from 
space. In th e space miss ions C0R0T (A uveranc et alJl2009h 
and Kepler (|Koch et al.ll20ld ). the photon noise determines 
the fundamental limit of photometric accuracy. However, 
the aperture of these instruments (C0R0T: 0.27 m, Kepler: 
0.95 m) are much smaller than the sizes of ground-based tele- 
scopes. This leads to a potential advantage of ground-based 
observations: it is possible to achieve a better temporal res- 
olution or to measure fainter objects. 

In the above example, the standard photometric B and 
R wide bands, often used for studies of the micro- variability, 
were adopted. Their large difference in wavelength helps to 
detect the colour variability of an astronomical object in the 
presence of spectral differences in light emitting processes. 
Asymptotic dependencies derived here allow us to obtain 
estimates for the opposite case of close and narrow photo- 
metric bands. For instance, to determine the variability in 
the spectral lines of radiation. When wavelengths difference 
is about 100 — 500 A, the effect of two decorrelation mecha- 
nisms (atmospheric dispersion and wind drift) is 1-2 orders 
of magnitude smaller than for the B and R bands. 
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6 CONCLUSION 

Although a large telescope aperture is comparable to the 
turbulence outer scale, we did not assess its potential im- 
pact on the obtained asymptotic dependencies. Probably, 
such accounting for can slightly change the dependence on 
the telescope diameter, but we do not expect a large effect 
because scintillation is produced mainly by the first and sub- 
sequent peaks of the Fresnel spectral filter. 

Also, implications of a non-Kolmogorov turbulence spa- 
tial spectrum were not ana lysed. The prelim inary estimate 
of this effect was made in (|Kornilovll20TTal 'l for the case of 
small apertures, where it was shown that a flatter spectrum 
weakens correlation of the scintillation in different photo- 
metric bands (increases colour scintillation). This is clear, 
because such spectrum leads to a raise of the contribution 
of the high frequencies which are transmitted by the differ- 
ential Fresnel filter. 

We summarise briefly the main results and conclusions 
obtained in this paper under the assumption of the Kol- 
mogorov spectrum of optical turbulence: 

- The colour scintillation was considered. It is produced 
by the wavelength dependence of the light diffraction in the 
turbulent atmosphere. 

- The spatio-chromatic scintillation caused by the atmo- 
spheric dispersion was analysed in addition. 

- The asymptotic dependencies were obtained to estimate 
the colour scintillation as a function of the telescope diam- 
eter. 

- The dependencies of the colour scintillation power on 
the turbulence altitude were determined (WFs for colour 
scintillation). 

- The colour scintillation is generated by highest turbu- 
lence and grows faster with the air mass than the normal 
scintillation. 

- Analysis of the behaviour of the scintillation was per- 
formed both for short exposures and for long exposures. The 
last is typical in conventional photometry. 

- In the short exposure regime, the achromatization of the 
scintillation occurs slowly. To halve the colour scintillation, 
the telescope diameter must increase by 8 times. 

- In the long exposure regime, the achromatization is hap- 
pening fast enough: when the telescope diameter increases 
by 2.3 times, the relative power of the colour scintillation is 
halved. 

- It was found that the combined effect of wind smoothing 
and atmospheric dispersion leads to a strong dependence of 
intensity of the colour scintillation on the wind direction in 
the upper atmosphere. 

- With an increase in the diameter of telescope, the colour 
scintillation is falling faster than the photon noise, so even- 
tually the limiting accuracy is determined by the photon 
noise. 
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Therefore, when observing on large telescopes, the scin- 
tillation in two different photometric bands is strongly cor- 
related. The colour scintillation is much less than the normal 
one and in the circumstances discussed above, the practical 
application of simultaneous multicolour measurements will 
provide distinct advantages in photometric accuracy. 



